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INTRODUCTION 


It  was  recognized  some  time  ago  that  vortex  interaction  effects  are  respon¬ 
sible  for  the  nonlinear,  high-lift  capability  of  slender,  low  aspect  ratio 
delta  wings.  As  a  consequence,  most  of  our  present  knowledge  of  vortex  inter¬ 
action  stems  from  extensive  studies  of  the  flowfields  of  such  planforms.  Re¬ 
cent  results,  though,  have  shown  that  wings  of  moderate  aspect  ratio  and  less 
sweepback  can  also  exhibit  high  lift  at  extreme  angles  of  attack  if  certain 
modifications  are  employed  to  create  favorable  vortex  interaction.  One  such 
modification  is  the  addition  of  a  canard  just  ahead  of  the  wing  in  the  so- 
called  close-coupled  canard  (CCC)  configuration. 

The  CCC  configuration  was  pioneered  in  Sweden  over  a  decade  ago  and  led  to 
the  successful  development  of  the  Saab  Viggen.  Interest  on  an  international 
scale  has  grown  since  then.  Recently,  Israel  demonstrated  a  new  version  of  its 
Kfir  fighter  which  incorporates  close-coupled  canards  to  improve  its  dogfight 
capability  in  air-to-air  combat.  Here  in  the  United  States,  the  Air  Force  is 
engaged  in  a  Control  Configured  Vehicles  Program,  using  canard  surfaces  on  a 
F-4E.  These  surfaces  are  part  of  conventional  and  unconventional  modes  of 
control  that  are  under  evaluation  to  enhance  the  effectiveness  of  this  air¬ 
craft.  The  Navy,  too,  has  devoted  considerable  attention  to  canards  for  future 
V/STOL  fighters,  like  the  XFV-12A . 

The  principal  advantages  that  the  CCC  configuration  brings  to  an  advanced  air- 
superiority  fighter  are  a  high  total  lift  capability  at  large  angles  of  attack, 
along  with  reduced  trim  drag.  These  aerodynamic  advantages  lead  to  impressive 
gains  in  combat  effectiveness  through  increased  turning  rates,  rate-of-climb 
and  acceleration.  Specific  excess  power,  from  an  energy  management  viewpoint, 
is  significantly  increased.  Other  attributes  include  improved  handling  qual¬ 
ities,  low  gust  sensitivities,  and  better  takeoff  and  landing  performance. 

The  existing  technology  has  been  sufficient  to  define  the  aerodynamic  potential 
of  the  CCC  configuration  for  a  range  of  specific  configurations.  A  fairly 
comprehensive  review  of  CCC  work  to  date  is  given  in  Ref.  1.  In  brief,  it  was 
almost  a  decade  after  the  studies  by  Saab  in  the  early  1960 's  before  further 
studies  were  pursued  in  the  United  States  at  NSRDC  and  NASA  Langley.  These 
studies  involved  mainly  wind  tunnel  tests  on  several  generalized  research 
models  and  a  modified  F-4B  model.  Prior  analytical  investigations  have  been 
very  limited,  usually  to  low  angles  of  attack  or  to  delta  planforms  for  both 
the  canard  and  the  wing.  Lately,  however,  NASA  Langley  has  been  trying  to  ex¬ 
tend  the  Polhamus  analogy  to  more  representative  planforms  at  high  angles  of 
attack . 

This  relatively  modest  effort  overall  has  not  provided  a  fundamental  under- 
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standing  of  the  flow  processes  that  contribute  to  the  superior  aerodynamic  per¬ 
formance.  That  is,  a  basic  description  of  the  favorable  vortex  interaction 
between  the  canard  and  the  wing  is  still  lacking,  especially  in  the  high  angle 
of  attack  regime.  Consequently,  it  is  not  surprising  that  there  is  no  analysis 
available  yet  to  predict  the  aerodynamic  performance  of  CCC  configurations. 

The  objective  of  the  present  investigation  is  to  obtain  a  good  physical  under¬ 
standing  of  the  canard/wing  flowfield  at  high  angles  of  attack  and  build  from 
this  a  simplified  theory  of  the  aerodynamic  interaction. 


APPROACH 


At  the  outset,  it  was  clear  from  existing  experimental  evidence  that  a  CCC  con¬ 
figuration  achieves  much  of  its  lift  at  high  incidence  from  vortex  lift.  This 
lift  is  generated  as  the  direct  result  of  flowfield  vortices  which  intimately 
interact  with  the  canard  and  wing  surfaces.  Knowledge  of  these  vortices  and 
corresponding  aerodynamic  data  are  essential  in  order  to  formulate  a  meaningful 
and  realistic  analysis.  Therefore,  a  comprehensive  investigation,  including 
experimental  and  theoretical  phases,  was  planned. 

The  four  guidelines  for  the  experimental  phase  of  the  investigation  were: 

■  Use  the  technology  developed  for  flow  visualization  with 
helium-filled  bubbles  to  identify  all  of  the  key  features 
of  the  interactive  flowfield. 

•  Select  a  particular  CCC  configuration  which  has  already 
demonstrated  a  significant  favorable  aerodynamic  inter¬ 
action. 

•  Make  the  model (s)  fairly  large  and  eliminate  or  minimize 
the  effects  of  any  fuselage  contribution. 

•  Obtain  a  consistent  set  of  force,  moment  and  pressure 
data  over  a  range  in  angle  of  attack  with  emphasis  on  high 
incidence. 

Two  separate  test  programs  were  necessary.  Both  were  conducted  in  the  NASA 
Langley  V/STOL  Research  Wind  Tunnel  which  has  a  test  section  14.5’  high,  21.75' 
wide  and  50 ‘  long.  It  can  be  run  either  closed,  slotted  or  open.  The  first 
test  program  concentrated  solely  on  visualization  of  the  flow  over  a  semi-span 
model,  mostly  with  the  test  section  open.  This  gave  an  ideal,  unobstructed 
view.  For  the  force,  moment  and  pressure  data  desired,  the  second  test  program 
utilized  a  sting-mounted,  full-span  model  with  different  wall  conditions.  The 
flow  visualization  results  and  wing  surface  pressure  data  constitute  important 
information  not  previously  known  for  any  CCC  configuration. 

For  the  theoretical  phase  of  the  investigation,  three  guidelines  were  esta¬ 
blished: 

•  Start  with  the  simplest  formulation  of  the  flowfield  that 
retains  the  key  features  observed  and  add  refinements  as 
required . 

•  Apply  potential  flow  theory,  incorporating  suitable 
singularity  distributions  of  vortex  lines  and  sheets 

for  the  canard  and  the  wing  surfaces  plus  any  free  vortex 
systems . 

•  Adopt  a  proven  numerical  method  as  far  as  possible  to 
solve  the  resulting  integral  equations. 


This  theoretical  approach  is  more  or  less  the  same  as  that  followed  historical- 


ly  to  attack  the  isolated  delta  wing  problem.  For  now,  the  method  must  be 
semi -empirical .  That  is,  the  strength  and  geometry  of  the  additional  free  vor¬ 
tex  systems  will  be  determined,  a  priori,  from  the  experiments.  Hopefully,  the 
insight  provided  will  eventually  lead  to  a  means  to  remove  the  need  for  empir¬ 
ical  data,  or  to  some  procedure  along  the  lines  of  the  Polhamus  analogy. 
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FLOW  VISUALIZATION 
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The  CCC  configuration  for  this  test  program  was  taken  from  one  of  the  NSRDC 
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generalized  research  configurations  .  It  is  representative  of  current  and  pro¬ 
jected  high-performance  aircraft  designs  and  exhibited  the  best  high-incidence 
lift  and  drag  characteristics. 

Fig.  1  shows  a  three-view  scale  drawing  of  the  half-span  model.  The  wing  has 
leading  and  trailing  edge  sweep  angles  of  50“  and  23.5°,  respectively,  with  a 
taper  ratio  of  0.151  and  an  aspect  ratio  of  3.90.  The  canard  is  half  of  a 
pure  60°  delta  planform,  giving  an  aspect  ratio  of  2.31.  During  most  of  the 
test,  the  canard  was  located  at  a  distance  equal  to  23%  of  the  wing  MAC  above 
the  wing  plane.  Longitudinally,  the  canard  trailing  edge  overlapped  the  wing 
by  21%  of  the  root  chord.  The  canard  surface  area  was  approximately  29%  of  the 
wing  area.  This  relatively  high  area  ratio  is  due,  primarily,  to  the  absence 
of  any  fuselage. 

The  wing  section  also  matched  the  NSRDC  wing  section,  a  NACA  64A008  airfoil 
laid  out  in  a  direction  inclined  25°  outward  relative  to  the  chordal  direction 
or  nearly  normal  to  the  trailing  edge.  For  the  canard,  a  flat  plate  section 
was  selected.  It  was  0.48%  thick  at  the  root  and  had  a  semi-circular  leading 
edge  and  blunt  trailing  edge. 

The  overall  size  of  the  model  was  dictated  by  flow  visualization  and  Reynolds 

number  considerations.  For  both,  the  largest  practical  size  was  desired.  The 

semi-span  finally  chosen  for  the  wing  was  6.0'  (1.83  m) .  This  gave  a  wing  MAC 
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of  43.66”  (1.11  m)  and  a  wing  reference  area  of  18.48  ft  (1.72  m  ).  Further 
description  of  the  model  is  contained  in  Ref.  1. 

The  use  of  helium-filled  bubbles  which  are  neutral ly-buoyant  have  proven  to  be 
an  accurate,  unique  tool  for  flow  visualization3.  Basically,  bubbles  about 
1/8”  (3  mm)  in  diameter  are  generated  by  a  small  probe,  or  head,  at  rates  up 
to  500  bubbles  per  second.  After  emerging  from  the  head,  the  bubbles  quickly 
reach  the  freestream  velocity  and  spread  out  substantially,  even  with  low 
levels  of  tunnel  turbulence.  If  the  flowfield  is  not  filled  adequately,  more 
bubbles  may  be  introduced  by  extra  heads.  Three  heads  were  adequate  in  this 
test  program. 

Still,  streak  photography  was  the  principal  means  for  recording  the  bubble 
trajectories.  With  a  high-sensitivity  vidicon,  videotape  recordings  were  also 
made  to  supplement  the  still  photography  and  identify  unsteady  effects.  In¬ 
tense  collimated  lighting  was  provided  by  a  searchlight  array  placed  downstream 
within  the  tunnel,  in  line  with  the  model. 
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Surface  flow  visualization  was  carried  out  in  order  to  examine  the  boundary 
layer  flows  on  the  wing  and  canard  as  well.  An  oil  flow  technique  with  tita¬ 
nium  dioxide  was  employed.  The  motions  of  the  oil  flow  patterns  were  recorded 
by  means  of  the  video  system. 

The  test  plan  was  to  concentrate  on  studying  the  model  flowfield  at  a  single 
angle  of  attack  in  the  regime  of  favorable  interaction.  From  the  force  data  in 
Ref.  2,  the  angle  decided  upon  was  a  =  25°.  To  compare  the  flow  pattern  at 
this  angle  with  that  at  a  much  lower  angle  where  the  wing  alone  would  not  be 
stalled,  a  =  10°  was  tested,  too.  Most  of  the  flow  visualization  was  done  at 
a  Reynolds  number  of  0.73  xlO6,  based  on  the  wing  MAC.  This  is  only  slightly 
less  than  the  value  for  the  NSRDC  data. 

Two  representative  streak  photographs  for  a  =  25°  are  presented  in  Fig.  2,  the 
upper  one  showing  the  wing  Without  the  canard  and  the  lower  one,  with  the  ca¬ 
nard.  Without  the  canard,  the  wing  is  badly  stalled  and  a  turbulent  separation 
region  covers  the  whole  upper  surface.  The  rest  of  the  flow  goes  smoothly 
around  the  separation  region,  which  is  fed  by  only  a  small  streamtube  close  to 
the  wing  root.  There  is  no  evidence  of  an  organized  leading-edge  vortex  or  a 
tip  vortex. 

Two  important  differences  in  the  wing  flowfield  are  apparent  when  the  canard  is 
installed.  First,  a  tight  leading-edge  vortex  is  formed  inboard  that  extends 
out  to  a  distance  about  equal  to  the  canard  semi-span.  This  vortex  is  not 
evident  in  Fig.  2  since  it  is  in  a  narrow  shadow  region.  Shown  dramatically, 
though,  is  the  bursting  of  this  vortex,  the  core  becoming  turbulent  and  en¬ 
larging  as  it  passes  over  the  wing  and  off  the  trailing  edge.  Outside  this 
turbulent  "funnel",  the  flow  is  laminar  but  it  has  appreciable  rotation. 

The  second  important  feature  is  the  unusual  ipanwiie  (£ou!  sandwiched  between 
the  burst  vortex  system  and  the  upper  surface  of  the  wing.  It  starts  inboard 
as  flow  which  has  gone  over  the  leading-edge  vortex  and  then  reattached  behind. 
After  traveling  underneath  the  burst  vortex,  it  almost  reaches  the  tip  before 
turning  back  again  into  the  freestream  direction.  No  tip  vortex  is  evident 
here  either. 

Although  not  shown  in  Fig.  2,  the  leading-edge  vortex  on  the  canard  bursts  near 
the  apex,  the  same  behavior  that  the  canard  alone  would  exhibit  at  this  angle 
of  attack.  The  core  of  this  burst  vortex  is  very  turbulent  and  grows  consid¬ 
erably  as  it  moves  rearward  to  the  trailing  edge  and  into  the  wake  above  Che 
wing.  It  appears  that  the  presence  of  the  wing  holds  this  burst  trailing  vor¬ 
tex  system  down  near  the  canard  and,  thus,  closer  to  the  wing  plane.  Between 
the  canard  and  wing,  the  flow  is  quite  smooth,  suggesting  a  "channeling”  effect 
which  preserves  the  leading-edge  vortex  on  the  inboard  portion  of  the  wing. 

The  surface  flow  visualization  distinctly  showed  the  separation  line  ahead  of 
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the  leading-edge  vortex  on  the  inboard  portion  of  the  wing.  In  many  other 
places,  however,  a  comparison  of  the  surface  patterns  with  the  outer  flow  pat¬ 
terns  revealed  significant  differences  in  flow  direction.  This  indicates  a 
drastic  change  in  flow  direction  through  the  boundary  layer. 

At  a  =  10®,  the  vortex  pattern  over  the  wing  without  the  canard  is  much  like 
the  pattern  at  a  =  25°  with  the  canard.  The  major  difference  lies  in  the  ab¬ 
sence  of  any  noticeable  spanwise  flow.  When  the  canard  is  added  at  this  low 
angle  of  attack,  there  are  no  pronounced  changes  in  the  wing  flowfield.  Just 
the  point  of  vortex  bursting  is  displaced  outboard  slightly. 

To  extract  an  average  or  quasi-steady  flowfield,  a  composite  picture  was  con¬ 
structed  by  superposition  of  streaks  from  a  number  of  photographs.  Consid- 
able  effort  was  needed  to  establish  the  three-dimensional  relationship  between 
the  streaks.  The  end  result  for  a  =  25®  is  drawn  in  Fig.  3.  The  broad  arrows 
describe  the  overall  flow,  which  is  steady  and  well  ordered.  The  narrow  arrows 
are  actual  bubble  trajectories  representative  of  the  random  motion  in  the  burst 
vortex  core. 

It  is  fairly  clear  that  the  flow  separates  all  along  the  wing  leading  edge, 
creating  a  leading-edge  vortex  sheet.  This  implies  that  there  is  no  pressure 
difference  or  load  across  the  wing  at  the  leading  edge,  equivalent  to  the 
Kutta-Joukowski  condition  at  the  trailing  edge.  Inboard,  the  leading-edge 
vortex  sheet  rolls  up  tightly  to  form  a  concentrated  vortex.  Outboard,  the 
sheet  wraps  around  the  turbulent  core  of  the  burst  vortex  as  the  core  seems  to 
detach  from  the  wing  surface.  This  causes  the  strong  spanwise  flow  that  was 
seen.  It  is  important  to  distinguish  this  high-velocity  spanwise  flow  outiidz 
the  boundary  layer  from  the  low-velocity  spanwise  flow  iniidz  the  boundary 
layer  on  yawed  wings  which  has  been  recognized  for  some  time.  As  a  consequence, 
lower  pressures  are  expected  oh  the  upper  wing  surface  outboard  and  so  more 
lift  from  the  tip  region  at  such  high  incidence. 

By  comparison,  the  inside  of  the  burst  vortex  is  nearly  a  dead-air  region.  The 
velocities  are  much  lower  and  almost  completely  random,  implying  that  the  pres¬ 
sure  is  practically  constant.  The  vortex  sheet  around  this  stagnant  region  is 
roughly  conical  in  shape  from  the  point  of  bursting  to  the  wing  trailing  edge. 
These  two  properties  of  the  burst  vortex  system  greatly  facilitate  the  theo¬ 
retical  formulation. 
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AERODYNAMIC  DATA 


The  full-span  model  for  the  second  test  program  was  built  to  the  same  scale  as 
the  semi-span  model,  or  12*  (3.66  m)  in  span.  Two  small  fuselages  housed  the 
internal  balances,  one  to  measure  the  total  forces  and  moments,  the  other  to 
measure  just  the  canard  loads.  The  upper  surface  of  the  right  wing  was  instru¬ 
mented  with  61  pressure  taps,  most  of  them  laid  out  at  7  spanwise  and  9  chord- 
wise  locations.  They  were  grouped  more  closely  near  the  leading  edge  in  the 
chordwise  direction. 

Full  six-component  force  and  moment  data  were  taken  by  means  of  two  sting 
knuckles  over  an  extended  range  in  a,  from  0°  to  36°.  Some  flow  visualization 
was  also  done  to  compare  with  the  previous  results.  Altogether,  four  config¬ 
urations  were  tested,  the  full  CCC  configuration,  the  wing-alone  or  wing  plus 
body,  the  canard-alone  or  canard  plus  body,  and  the  body  by  itself.  In  addi¬ 
tion,  the  tunnel  velocity  was  varied  to  assess  the  effects  of  Reynolds  number. 
The  Reynolds  number  based  on  the  wing  MAC  ranged  from  approximately  0.7  x  106 
to  3.7  x 106. 

The  lift  curve  for  the  full  CCC  configuration  is  plotted  in  Fig.  4.  It  is  seen 
that  the  lift  coefficient  increases  monotonically ,  though  not  linearly,  up  to 
a  value  of  over  1.8  referenced  to  the  wing  area.  There  is  no  indication  of 
stall  up  to  this  point. 

To  evaluate  the  favorable  aerodynamic  interaction  between  the  canard  and  wing , 
a  second  set  of  data  points  has  been  plotted  for  comparison.  These  points  are 
for  an  "uncoupled"  canard  configuration  in  which  the  canard  and  wing  act  in¬ 
dependently  of  each  other.  They  were  obtained  by  adding  the  lift  coefficients 
for  the  wing-alone  and  canard-alone  configurations,  then  subtracting  the  body- 
alone  data.  Only  data  for  high  angles  of  attack  are  shown  because  canard-alone 
and  body-alone  data  were  not  taken  at  low  angles.  This  uncoupled  canard  con¬ 
figuration  stalls  out  between  25°  and  30° . 

As  found  by  others,  the  curves  cross  over.  Favorable  interaction  takes  place 
above  20°,  the  lift  increment  between  the  close-coupled  and  the  uncoupled 
canard  configuration  increasing  with  angle  of  attack.  Below  20°,  the  uncoupled 
lift  is  a  little  higher. 

This  behavior  can  be  explained  by  examination  of  the  individual  contributions 
of  the  wing  and  canard.  The  wing-alone  data  are  shown  in  Fig.  5.  The  initial 
lift  curve  slope  for  the  wing  in  the  CCC  configuration  is  appreciably  less  than 
for  the  wing  alone.  This  should  be  anticipated  from  classical  wing  theory  be¬ 
cause  of  the  downwash  at  the  wing  from  the  canard.  By  16°,  the  lift  curve 
slopes  are  about  equal.  The  wing-alone  lift  peaks  out  at  25°,  while  the  lift 
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Figure  4.  Comparison  Of  Close-Coupled  Canard  And 
Uncoupled  Canard  Lift  Data 


+  WING  IN  CCC  COFIGURATION  -  LOW  ALPHA 


generated  by  the  wing  in  the  presence  of  the  canard  continues  to  cirmb,  al¬ 
though  at  a  declining  rate.  In  boundary  layer  terminology,  the  canard  appears 
to  inhibit  or  delay  stalling  of  the  wing  at  high  angles  as  substantiated  by 
Fig.  2. 

The  canard-alone  data  are  shown  in  Fig.  6,  the  lift  coefficients  again  refer¬ 
enced  to  the  wing  area.  They  can  be  multiplied  by  the  wing/canard  area  ratio 
of  3.399  to  refer  them  to  the  canard  area  if  desired.  Extrapolating  these  data 
to  the  origin,  the  lift  curve  slope  there  is  less  than  for  the  canard  in  the 
CCC  configuration.  As  before,  this  should  be  anticipated  from  classical  wing 
theory  because  of  the  upwash  at  the  canard  from  the  wing.  Over  the  entire 
angle  of  attack  range,  the  lift  and  lift  curve  slope  for  the  canard  in  the 
presence  of  the  wing  is  greater  than  for  the  canard  alone.  At  30°,  the  canard- 
alone  lift  reaches  a  maximum,  indicating  stall,  but  the  canard  lift  in  the  CCC 
configuration  still  increases.  Clearly,  the  improvement  in  canard  lift  at  low 
angles  from  favorable  interaction  with  the  wing  does  not  fully  offset  the  re¬ 
duction  in  wing  lift.  At  high  angles,  though,  the  interaction  is  favorable 
to  both  surfaces. 

It  is  significant  to  note  that  the  canard  develops  more  lift  above  20°  coupled 
to  the  wing  than  the  maximum  lift  which  the  canard  alone  develops.  There  is 
no  way  to  account  for  this  by  a  change  in  the  effective  angle  of  attack  of  the 
canard  due  to  the  upwash  from  the  wing.  Instead,  a  change  in  the  effective 
camber  would  have  to  be  introduced  or,  perhaps,  something  more  profound. 

It  is  also  significant  to  note  that  the  high-angle  data  of  Fig.  6  for  the  ca¬ 
nard  in  the  CCC  configuration  can  be  fitted  extremely  well  by  the  Polhamus  lift 

4 

equation  for  delta  wings  .  If  the  coefficients  are  adjusted  to  match  the  mea¬ 
sured  values  of  CL  at  16°  and  32°,  the  equation 

2  2 

C^  =  0.974  sin  a  cos  a  +  0.777  sin  a  cos  a  (1) 

predicts  the  measured  values  at  all  other  angles  of  attack  from  14°  to  34° 
within  0.5*.  Multiplying  the  constants  in  Eq.  (1)  by  3.399,  they  can  be  com¬ 
pared  with  the  values  of  Kp  and  Kv  in  Ref.  5  for  a  60°  delta  wing.  The  poten¬ 
tial  lift,  or  first  term  above,  is  37%  higher  than  given  by  the  Polhamus 
analogy  for  an  isolated  delta  wing.  On  the  other  hand,  the  vortex  lift,  or 
second  term  above,  is  18*  lower.  Nonetheless,  the  fact  that  this  lift  equation 
derived  from  the  Polhamus  analogy  is  so  good  is  very  encouraging  from  a  theo¬ 
retical  viewpoint. 

With  respect  to  Reynolds  number  effects,  the  lift  data  were  repeatable  over  the 
whole  range  in  Reynolds  numbers  covered  in  the  test  program.  This  is  also  en¬ 
couraging  because  it  indicates  that  any  flowfield  model  which  is  adopted  can  be 
used  with  some  generality. 
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The  flow  visualization  results  agreed  on  the  whole  with  the  results  from  the 
first  test  program.  Only  the  geometry  of  the  burst  vortex  on  the  wing  was 
slightly  different.  The  included  angle  was  somewhat  larger.  Also,  the  up¬ 
stream  edge  of  the  burst  vortex  moved  forward  and  became  almost  parallel  to  the 
leading  edge.  At  the  two  angles  of  observation,  25®  and  32°,  the  leading  edge 
vortices  on  the  canard  were  fully  burst  as  seen  in  the  semi-span  tests.  This 
may  be  the  reason  for  the  lower  vortex  lift  from  Eq.  (1). 

A  sample  of  the  pressure  data  is  illustrated  in  Fig.  7  as  contour  lines  of  con¬ 
stant  pressure  coefficient  over  the  upper  surface  of  the  wing.  Contours  for 
coefficient  values  of  -0.5,  -1,  -2,  -3  and  -4  were  chosen.  Inboard,  the  con¬ 
tours  are  close  together  near  the  leading  edge  under  the  leading  edge  vortex. 
Outboard,  the  contour  lines  for  -1  and  -2  are  "pulled"  outward  and  rearward  by 
the  burst  vortex  system.  The  extreme  points  of  these  curves  lie  almost  di¬ 
rectly  under  the  centerline  of  the  burst  system.  These  low  pressures  must  be 
created  by  the  high-velocity  spanwise  flow.  Inboard  or  outboard,  the  shape  of 
the  contour  lines  do  not  change  a  great  deal  over  the  range  of  high  angles  of 
attack.  This  shows  once  more  the  persistence  of  the  flowfield  pattern. 

Some  comments  should  be  made,  finally,  regarding  tunnel  wall  effects.  At  the 
suggestion  of  NASA,  the  full  CCC  configuration  was  tested  with  the  wind  tunnel 
test  section  closed,  slotted  and  open.  The  closed  test  section,  which  was 
used  in  all  other  runs,  gave  the  highest  lift  values.  With  the  slotted  test 
section,  both  the  total  lift  and  the  canard  lift  were  4%  to  6*  lower  over  the 
high-angle  range.  With  the  open  test  section,  i.e.  raising  the  two  walls  and 
ceiling,  these  lift  values  were  9*  to  12%  lower.  Wall  corrections  are  evi¬ 
dently  needed,  but  the  normal  corrections  did  not  give  reasonable  answers.  The 
theoretical  analysis  should  enable  the  derivation  of  appropriate  corrections 
for  this  configuration. 

The  canard  lift  values  with  the  tunnel  slotted  and  open  were  also  fitted  very 
well  by  the  Polhamus  lift  equation.  Surprisingly,  the  potential  lift  of 
Eq.  (1)  for  the  test  section  closed  was  essentially  unchanged.  The  lift  re¬ 
duction  in  both  cases,  therefore,  came  entirely  in  the  vortex  lift.  Compared 
to  the  theoretical  value  for  an  isolated  delta  wing,  as  before,  the  vortex 
lift  dropped  from  a  value  18%  lower  with  the  closed  test  section  to  35%  lower 
with  the  slotted  test  section  and  down  to  47%  lower  with  the  open  test  sec¬ 
tion. 


16 


CCC  THEORY 


To  construct  a  simplified  theory  for  such  a  complex  and  unusual  flowfield  is  a 
formidable  challenge.  Before  proceeding,  though,  it  is  necessary  to  examine 
two  aspects  of  the  flowfield  first  in  some  detail.  One  is  the  spanwise  flow 
and  the  other,  the  burst  vortex  system  above  the  wing. 

The  upper  sketch  of  Fig.  8  illustrates  the  velocity  vectors  for  a  representa¬ 
tive  point  on  the  wing  and  the  associated  vorticity  required  to  replace  the 
wing  by  a  vortex  sheet.  This  point  is  located  at  the  mid-chord,  directly  under 
the  centerline  of  the  burst  vortex.  The  vector  lies  in  the  direction  of  the 
spanwise  flow  above  the  wing,  while  g2  lies  in  the  direction  of  the  chordwise 
flow  observed  directly  under  the  wing.  The  relative  magnitude  of  these  vectors 
were  determined  from  the  upper  surface  pressure  coefficient  of  Fig.  7  and  an 
estimate  of  the  lower  surface  pressure  coefficient.  The  "sheet"  velocity 
is  the  mean  of  g.  and  g2 . 

The  local  vorticity  required  to  induce  this  change  in  velocity  across  the 
sheet  is  given  by  the  vector  equation 

1  =  -  (gx  -g2>  x  n  (2) 

where  n  is  the  upward  facing  unit  vector  normal  to  the  wing  plane.  Surprising¬ 
ly,  ^  is  pointed  almost  in  the  chordwise  or  streamwise  direction.  The  bound 
vorticity  in  classical  wing  theory  is  usually  pointed  in  the  spanwise  direc¬ 
tion. 

In  a  similar  fashion,  the  jump  conditions  can  be  found  for  a  vortex  sheet 
around  the  burst  vortex  that  will  insure  essentially  stagnant  flow  inside.  The 
vectors  involved  in  this  case  are  shown  in  the  lower  sketch  of  Fig.  8  for  a 
point  on  the  upper  side  of  the  burst  system  taken  above  the  previous  point. 

The  velocity  vector  g^  outside  the  sheet  was  determined  by  flow  visualization. 
If  =  0  immediately  inside  the  sheet,  gg  =  g^/2  and  follows  from  Eq.  (2)  . 
The  required  vorticity  here  on  the  sheet  has  components  both  in  the  circum¬ 
ferential  and  ray  directions  on  the  burst  vortex. 

With  this  understanding  in  hand,  a  theoretical  model  can  be  formulated.  This 
model  is  drawn  in  schematic  form  in  Fig.  9.  For  the  canard  and  wing,  conven¬ 
tional  horseshoe  vortices  are  used.  The  burst  vortex  system  above  the  wing  is 
idealized  as  a  cone,  composed  of  ring  vortices  plus  "cranked"  horseshoe  vor¬ 
tices.  The  ring  vortices  account  for  the  circumferential  component  of  the 
vorticity  vector  on  the  burst  vortex  sheet.  The  cranked  horseshoes,  depicted 
in  Fig.  9  by  a  single  horseshoe,  are  composed  of  three  parts.  On  the  burst 
vortex  sheet,  the  conical  elements  account  for  the  ray  component  of  the  vor- 


ticity  vector.  These  elements  are  smeared  around  the  cone  but  are  fed  at  only 
one  azimuth  by  the  umbilical  elements  from  the  wing  leading  edge.  The  umbil¬ 
ical  elements,  in  turn,  join  the  wing  plane  elements  which  may  be  skewed  some¬ 
what  inboard  or  outboard.  The  wing  plane  elements  are  very  important  to  vortex 
lift,  as  discussed  later  on. 

The  potential  flow  for  present  purposes  will  be  defined  as  the  regular  attached 
flow  in  the  absence  of  the  burst  vortex  system.  To  determine  this  flow,  the 
normal  component  of  velocity  must  vanish  everywhere  on  the  wing  and  the  canard. 
If  the  surfaces  are  parallel  and  uncambered,  it  follows  then  that 

U  sin  “  +  www  +  wwc  =  0  (3) 

U  sin  a  +  wcc  +  wcw  =  0  (4) 

The  velocity  w^  is  the  upwash  on  the  wing  induced  by  the  vortex  system  of  the 

wing,  wwc  is  the  upwash  on  the  wing  induced  by  the  vortex  system  of  the  canard, 

wcc  is  the  upwash  on  the  canard  induced  by  the  vortex  system  of  the  canard  and 
wcw  is  the  upwash  on  the  canard  induced  by  the  vortex  system  of  the  wing. 

Eqs.  (3)  and  (4)  are  exact  for  arbitrary  a  within  this  formulation. 

These  equations  constitute  a  set  of  integral  equations  coupled  through  wwc  and 
wcw>  They  can  be  effectively  decoupled  and  solved  in  closed  form  by  successive 
iteration  using  operators.  The  final  solutions  may  be  expressed  as 


0  v  {y  (0)  +  y  (1)} 

w  Tw  Tw  ; 


Yc  =  l  °cV  Wc<0)  +  Yc<1>} 
°  V=0  C  c  c 


where  yw  is  the  strength  of  the  unknown  horseshoe  vortex  distribution  on  the 
wing  and  yc>  the  strength  on  the  canard.  The  operators  0  and  0c  are  fairly 
complex  operators  given  by  the  product  of  four  other  operators ,  two  corre¬ 
sponding  to  the  inversion  of  Eqs.  (3)  and  (4)  with  no  interaction  and  two,  to 
the  upwash  of  these  surfaces  on  each  other.  With  respect  to  the  remaining 

terms,  y  ^  is  the  vortex  strength  on  the  isolated  wing  and  y  ,  on  the 

w  (1)  (1)  c 
isolated  canard.  The  other  vortex  strengths,  yw  '  and  yc  ,  arise  from  just 

the  upwash  induced  by  the  vortex  distributions  for  the  isolated  surfaces. 

Any  one  of  several  computational  methods  for  lifting  surfaces  in  potential  flow 
could  be  utilized  to  carry  out  the  solutions  of  Eqs.  (5)  and  (6)  numerically. 
The  Douglas  EVD  Method®'^  was  selected  here.  Since  it  was  limited  to  vortex 
panels  in  a  single  plane,  this  method  had  to  be  generalized  to  include  the 


influence  of  out-of -plane  panels. 


The  calculations  for  the  influence  of  the  burst  vortex  sheet  are  considerably 
simplified  if  the  induced  flowfield  is  conical.  This  can  be  achieved  by  as¬ 
suming  a  constant  distribution  of  vorticity  and  extending  the  conical  sheet 
from  the  apex  to  infinity.  Then  the  induced  velocities  are  constant  along  rays 
or  concentric  cones . 

A  cylindrical  coordinate  system  (x,r,0)  is  employed  with  the  origin  at  the  apex 
and  the  x-axis  along  the  cone  axis.  The  circumferential  vorticity  y0 ,  or  dis¬ 
tribution  of  ring  vortices,  induces  an  axial  velocity  component  u  and  a  radial 
velocity  component  v.  The  0-integration  of  the  Biot-Savart  law  can  be  per¬ 
formed  analytically  in  terms  of  Legendre  functions  of  the  second  kind  and  half¬ 
integer  order  Qn_^/2'  the  conical  variables  n  =  x/r  and  nv  =  xv/rv  are  usec* 

to  indicate  the  rays  on  which  the  field  point  (x,r)  and  vortex  points  (xv,ry) 
are  located,  the  results  can  be  reduced  to  the  form  shown  in  Table  1.  On  the 
other  hand,  the  ray  vorticity  y^  induces  only  a  circumferential  velocity  com¬ 
ponent  w.  It  can  be  reduced  like  u  and  v  and  is  included  along  with  them  in 
Table  1. 

This  leaves  the  umbilical  elements  and  the  wing  plane  elements  of  the  cranked 
horseshoe  vortices.  Since  the  influence  of  the  umbilical  elements  should  be 
small,  these  elements  will  be  neglected.  The  wing  plane  elements,  though, 
will  play  a  pivotal  role.  This  role  will  be  discussed  further  in  a  moment. 

With  regard  now  to  numerical  results  for  the  potential  lift,  calculations  were 
performed  for  the  canard  alone,  wing  alone  and  the  full  CCC  configuration. 

These  were  run  on  a  Hodcomp  11/220  Computer,  taking  a  Glauert  or  cosine  spacing 
of  the  EVD  panels  in  both  the  spanwise  and  chordwise  directions.  In  general, 
the  agreement  with  other  theoretical  and  experimental  data,  as  well  as  with  the 
present  experimental  data,  is  very  good. 

For  the  canard  alone,  the  computed  values  of  CL  /sin  a  range  from  2.51  to  2.62. 
These  values  correspond  to  12  x 6  and  6x6  arrays  of  EVD  panels,  respectively, 
where  the  first  integer  is  the  number  of  panels  counted  in  the  spanwise  di¬ 
rection  and  the  second  integer,  in  the  chordwise  direction.  Values  from 
Refs.  4  and  5  for  other  lifting  surface  methods  are  2.45  and  2.41,  by  compari¬ 
son,  or  about  5%  less.  The  experimental  values  reported  in  Ref.  8  for  low 
angles  are  2.58  and  2.61,  or  exceptionally  close  to  the  theoretical  values, 
while  the  value  of  2.75  from  Ref.  9  is  a  little  high. 

For  the  wing  alone,  the  computed  values  of  CL  / sin  a  range  from  3.29  to  3.50. 
The  lowest  value  corresponds  to  an  18  x 6  array  and  the  highest  value,  to  a  6x6 
array.  The  present  experimental  data  at  low  angles,  see  Fig.  5,  fall  between 
3.44  and  3.55  with  the  test  section  closed.  This  is  somewhat  higher  than  the 
theoretical  predictions,  as  might  be  expected  from  the  tunnel  interference 
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discussed  earlier.  An  estimate  of  the  comparable  values  for  a  slotted  test 
section  are  3.28  to  3.39,  agreeing  quite  well  with  the  calculated  results.  A 
value  of  3.39  was  obtained  from  the  NSRDC  data  in  Ref.  2. 

Computer  runs  for  the  full  CCC  configuration  were  made  with  two  sets  of  EVD 
panel  layouts.  One,  with  a  12  x 6  array  on  the  wing  and  a  6x6  array  on  the  ca¬ 
nard,  gave  a  value  for  /sin  a  based  on  the  wing  area  of  3.71.  The  other, 
with  15x6  and  9x6  arrays,  instead,  gave  a  slightly  higher  value  of  3.73.  By 
comparison,  the  present  test  data  at  low  angles  varied  from  3.80  to  3.95  with 
the  test  section  closed  and  from  3.64  to  3.73  with  the  test  section  slotted. 

Numerical  results  for  the  velocities  induced  by  the  burst  vortex  sheet  u/y8# 
v/y0  and  w/y^  are  plotted  in  Fig.  10.  These  values  were  found  by  numerical 
integration  of  Eqs.  (7),  (8)  and  (9)  for  a  semi-vertex  angle  of  11s  as  measured 
from  flow  visualization  photographs.  Careful  examination  of  Eq.  (7)  reveals 
that  this  integral  is  actually  slowly  divergent.  Therefore,  the  numerical  in¬ 
tegration  was  cut  off  at  £  *=  50. 

The  position  of  the  vortex  sheet  ny  is  denoted  by  the  dashed  line  in  Fig.  10. 
Directly  outside  the  sheet,  the  dimensionless  axial,  radial  and  circumferential 
velocities  are  -0.128,  0.191  and  1.000,  respectively.  Inside  the  sheet,  u/y0 
is  approximately  equal  to  -1.110  everywhere  and  the  other  two  components 
identically  vanish.  These  values  satisfy  the  correct  jump  conditions  across 
the  sheet  of  cos  11®,  sin  11®  and  1.000,  respectively. 

To  qualify  the  theoretical  model  of  the  burst  vortex  system  up  to  this  point, 

particle  trajectories  were  computed  for  the  flow  about  the  conical  vortex  sheet 

placed  in  a  uniform  stream.  The  magnitude  of  the  stream  velocity  was  chosen 

to  cancel  u  inside  the  sheet  and  the  ratio  of  y„  to  y_  was  taken  from  the  flow 

o  n 

inclination  measured  just  outside  the  upper  surface  of  the  burst  vortex,  see 
Fig.  8.  Under  these  conditions,  the  conical  vortex  sheet  simulates  a  solid 
cone  so  there  is  no  flow  through  the  sheet. 

Several  different  "starting  points"  were  tried  for  the  "particles"  and  the  tra¬ 
jectories  superimposed  upon  photographs  of  the  actual  bubble  traces.  In  cer¬ 
tain  areas,  some  local  discrepancies  exist,  such  as  close  to  the  upper  wing 
surface.  Overall,  though,  the  shape  and  behavior  in  the  two  cases  are  remark¬ 
ably  alike. 

Once  the  flowfield  model  for  the  burst  vortex  system  has  been  established,  the 
final  steps  of  development  of  the  theory  can  be  undertaken.  The  first  thing 
that  has  to  be  done  is  the  computation  of  the  additional  downwash  contribu¬ 
tions.  With  empirical  values  for  ya  and  y_»  the  contributions  of  the  ring 

o  n 

vortices  and  the  ray  elements  of  the  cranked  horseshoe  vortices  can  be  readily 
calculated  from  u,  v  and  w.  Since  the  contribution  of  the  umbilical  elements 
of  the  cranked  horseshoe  vortices  will  be  neglected,  only  the  contribution  of 


the  wing  plane  elements  is  needed.  This  can  be  carried  out  in  the  manner  of 
lifting-line  theory  after  the  skew  angle  has  been  fixed. 


* 

With  the  additional  downwash  contributions  computed,  w^  and  wcw  in  Eqs.  (3) 
and  (4)  can  be  adjusted  and  new  distributions  of  and  yc  determined  from 
Eqs.  (5)  and  (6).  The  operators  0v  and  0c  will  have  to  be  changed  for  this 
purpose  by  suppressing  the  leading  edge  singularity.  Also,  the  numerical  im¬ 
plementation  will  have  to  be  changed,  but  the  explicit  identification  of  this 
singularity  in  the  Douglas  EVD  Method  should  facilitate  this  change. 

The  last  step  of  development  of  the  theory  will  be  the  computation  of  the  sur¬ 
face  pressures  and  the  resultant  lift.  To  do  this,  the  total  surface  veloc¬ 
ities  will  be  found,  then  substituted  into  the  full  Bernoulli  equation  and 
integrated.  It  can  be  seen  immediately  that  the  "lift  mechanism*  will  be  quite 
different.  That  is,  besides  the  lift  produced  by  the  product  of  the  X-compo- 
nent  of  the  surface  sheet  velocity  and  the  Y-component  of  the  vorticity,  there 
will  be  lift  due  to  the  product  of  the  Y-component  of  the  surface  sheet  ve¬ 
locity  and  the  X-component  of  vorticity. 

This  latter  lift  depends  strongly  upon  the  wing  plane  elements  of  the  cranked 
horseshoes  vortices.  In  fact,  it  constitutes  most  of  the  vortex  lift,  or  lift 
increment  in  the  presence  of  the  burst  vortex  system.  Judged  from  Fig.  8,  the 
vortex  lift  will  dominate  the  potential  lift  over  much  of  the  wing  at  high 


CONCLUSIONS 


Based  upon  the  progress  to  date  of  a  comprehensive  investigation  to  develop  a 
simplified  theory  for  the  close-coupled  canard  configuration,  the  following 
conclusions  are  drawn: 

•  A  detailed  understanding  of  the  flowfield  for  a  represent¬ 
ative  configuration  has  been  obtained  on  a  6'  (1.83  m)  semi¬ 
span  model  by  flow  visualization. 

•  Three  prominent  features  of  the  flow  have  been  identified. 

One  is  the  presence  of  a  large,  well-ordered  burst  vortex 
above  the  wing  at  high  angles  of  attack.  Another  is  the 
strong  spanwise  flow  which  is  sandwiched  between  this 
vortex  and  the  wing  at  these  angles.  The  third  is  the  ap¬ 
parent  constraint  of  the  canard  vortex  system  that  greatly 
improves  the  canard  lift  characteristics. 

•  Corresponding  aerodynamic  measurements  have  been  taken  on 
a  12'  (3.66  m)  full-span  model,  including  pressure  con¬ 
tours  on  the  upper  surface  of  the  wing .  These  contours 
clearly  correlate  with  the  burst  vortex  observed. 

•  Highly-favorable  interaction  was  achieved  with  lift  values 
over  1.8  reached  at  a  =  34°  without  stall.  The  canard 
lift  has  been  fitted  within  0.5*  variation,  over  the  whole 
high-angle  range  tested,  with  the  Polhamus  lift  equation 
by  adjusting  the  coefficients. 

•  The  proper  velocity/vorticity  relationships  for  both  the 
wing  and  the  burst  vortex  have  been  derived  and  a  math¬ 
ematical  model  has  been  formulated.  This  involves  besides 
the  usual  singularities  the  addition  of  ring  vortices  and 
cranked  horseshoe  vortices  to  represent  the  burst  vortex 
system. 

•  Two  other  key  steps  have  been  completed.  The  potential 
lift  has  been  found  in  closed  form  for  arbitrary  angles  of 
attack.  Also,  the  influence  of  the  most  difficult  parts 
of  the  burst  vortex  system  has  been  reduced  to  a  conical 
form  which  is  very  convenient  for  computations. 

•  Both  steps  have  been  qualified  numerically.  Lift  results 
calculated  for  low  angles  by  extending  the  Douglas  EVD 
Method  show  excellent  agreement  overall  with  other  theo¬ 
retical  and  experimental  data,  as  well  as  with  the  present 
experimental  data.  Particle  trajectories  of  the  flow 
about  the  idealized  burst  vortex  in  a  uniform  stream  have 
been  computed,  too,  and  compare  well  with  the  actual  tra¬ 
jectories  seen. 

•  At  high  angles  of  attack,  the  vortex  lift  clearly  dominates 
the  potential  lift.  Some  work  remains  to  be  done,  but  the 
mechanism  for  this  lift  has  been  determined. 
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1  3  ABSTR AC  T 


The  close-coupled  canard  configuration  has  emerged  as  an  advanced 
design  for  highly-maneuverable  fighter  aircraft,  achieving  much  of 
its  lift  at  high  angles  of  attack  from  vortex  lift.  This  paper  sum¬ 
marizes  the  progress  to  date  on  a  comprehensive  investigation  to 
develop  a  simplified  theory  for  the  aerodynamic  performance.  Flow 
visualization  tests  on  a  6'-- (1.83  m)  semi-span  model  of  a  representa¬ 
tive  configuration  are  first  discussed.  These  tests  are  followed  by 
measurements  on  a  12'  (3.66  m)  full-span  model  to  obtain  force  and 
moment  data,  as  well  as  pressure  data  on  the  upper  surface  of  the 
wing.  With  this  understanding,  the  proper  velocity /vorticity  rela¬ 
tionships  on  the  wing  and  the  burst  vortex  above  the  wing  are  found 
and  a  compatible  mathematical  model  is  formulated.  The  general 
solution  to  the  potential  flow  problem  is  then  given  and  the  influ¬ 
ence  of  conical  parts  of  the  burst  vortex  system  is  derived.  Numeri¬ 
cal  results  for  these  steps  look  very  good.  In  conclusion,  the  final 
steps  in  the  development  of  the  theory  are  outlined  and  a  new  mech¬ 


anism  for  the  vortex  lift  is  identified. 
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